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^y \ Abstract. With a generally covariant equation of Dirac fields outside a black 

hole, we develop a scattering theory for massive Dirac fields. The existence of 
^3 ■ modified wave operators at infinity is shown by implementing a time-dependent 

Q \ logarithmic phase shift from the free dynamics to offset a long-range mass term. 

O^' The phase shift we obtain is a matrix operator due to the existence of both 



^H . positive and negative energy wave components. 
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1. Introduction 

As a classical scattering theory of massive Dirac fields on black hole, this paper 
is devoted to proving the existence of modified wave operators with logarithmic 
phase shifts at infinity, along the line of the treatment of the Coulomb problem 
by Dollard [1]. Generally speaking, the classical scattering off a black hole, on its 
own, does not have much physical meaning since it is essentially not observable. 
However, it is a first and necessary step towards the better understanding of 
relativistic quantum Dirac fields on black hole, which is investigated by the author 
elsewhere [2]. Mathematically, it is interesting since it involves some powerful 
mathematical techniques. 

The time-dependent scattering theory for scalar fields on the Schwarzschild 
metric has been established by Dimock and Kay [3] [4] [5] [6]. Some recent 
investigations of asymptotic completeness and existence have been taken mainly 
by Bachelot on the Klein-Gordon equation [7] and on the Maxwell system [8], 
and by Bachelot and Motet-Bachelot on the black- hole resonance [9]. 

For Dirac fields, there are some early discussions on the asymptotic behavior 
of Coulomb scattering by Dollard and Velo [10], and by Enss and Thaller [11]. 
On black hole, the existence of classical wave operators for Dirac fields has been 
proved at horizon (with and without mass) and at infinity (without mass) by 
Nicolas [12]. We proceed to clarify the asymptotic problem at infinity for massive 
Dirac fields, and show that the long-range mass term in the Hamiltonian can 
actually be offset by a logarithmic phase shift from the free dynamics along the 
line of the work for massive scalar fields by Dimock and Kay [3]. 

This paper is organized as follows. In next section, we start by showing 
some standard results, especially the relativistic Dirac equations on black hole, 
which are essential to the development of our theory. After acquiring enough 
preliminary knowledge, we develop this scattering theory of massive Dirac fields 
on black hole at infinity in two steps. As a first step, to get a sense of physical 
process, we give a rough qualitative picture by considering the propagation of a 
peaked wave in Section 3. Then, as a second step, we use some mathematical 
tools to prove the existence of modified wave operators for any smooth wave in 
Section 4. Finally in Section 5 we add some remarks that are helpful to our 
understanding and future research. 

2. Preliminaries 

In the presence of gravitational field, the space-time structure is considered as 
pseudo-Riemannian. However we can still fix the locally inertial frame in which 



the spin representation of Lorentz group can be established. The vierbein fields 
are any collection of four vectors V^{x) which satisfy 

< Va, Vb >= g^,,{x)V^^{x)Vi'{x) = rjab, (1) 

with a Minkowskian metric j ??a6 r = diag{l, —1, —1, —I). Here Greek indices 
denote space-time components and Latin indices represent spinor components. 
Clearly (1) is equivalent to 

g^^ix) = V,^ixW{x)rj'^\ (2) 

The relation 

g^''{x) = \{j>^ix),Yix)}, (3) 

can be satisfied by choosing spinor-tensors 

7^(x) = V,^{x)j\ (4) 

where 7" are Dirac matrices with the convention in [13]. The vierbein fields 
V^{x) transform like coordinate vectors and local Lorentz vectors 

V,^{x)=Li{x)V,^{x), (5) 

with a map L : M ^ £(Lorentz group). Then Dirac spinor fields ip{x) : M -^ C^, 
transform as 

V'(x) ^ S{L{x))ij{x), (6) 

where S is the finite spinor representation of the Lorentz group. 

The next goal is to define generally covariant derivatives that render the 
transformation (6) invariant. This is achieved by the requirement that j'^ix) 
have vanishing covariant derivatives [14] [15] 

V;.7'^ = I?^7"-r^7" + 7% = 0. (7) 

Then spin affine connections can be determined as [16] [2] 

T, = ^G^-''\V,V:)V,,, (8) 

where V^ are the covariant derivatives in space-time only, and the generators of 
the spin representation of Lorentz group are 

GM=l[^a,^fe]. (9) 



With the introduction of spin affine connections, we have the covariant derivatives 
V^ = d^ — Tfi on Dirac spinor fields [2]. Finally, we use the Equivalence Prin- 
ciple of Relativity to generalize the Dirac equation by simply replacing regular 
derivatives with covariant derivatives: 

ijf'idf, - r^)^ = m^. (10) 

Now let us consider a spherical black hole with a metric 

{g^,}=diag{(,-r\-r^-{rsmef), (11) 

where .^ = 1 — 2GM/r. From (1), the vierbein fields on such a black hole can be 
chosen as 

{V,^} = d^agir'^^e^^r-\irsme)-'). (12) 

Then the spin affine connection matrices of a black hole are [15] [2] 

r, = i\e^VA le^Vl', le^'sine^h' + \cos9^V). (13) 

Now making a substitution in spherical coordinates, 

^ = e^H-g)'^'-^, (14) 

where g = det{g'^'^) = — r^(sin^)^, one has a nice equation 

i^^d^^ = m^. (15) 

Rewrite it as follows [14] [2] 

idt^ = H^, H = ai® [ai{-i^—) + ^K] + (ag /)me^/^ (16) 

or r 

where Cj are Pauli matrices, / is a 2 x 2 unit matrix, represents a tensor 
product, and K is an angular operator 

.^ • ^ • 1 ^ 

-^ = ~^^2 757-*o-3— - — . (17) 

oO sm ocp 

By using the Regge- Wheeler variable r^ defined by dr^/dr = ^~^, the Hamiltonian 
in (16) can also be written as 

d fi/2 

Hr, = ai [ai{-i^—) + ^ — K] + (cjg I)mE}/^. (18) 

or^ r 

Later in Theorem 1, we will prove the self-adjointness of this Hamiltonian. 



According to Stone's theorem [17], the self-adjointness of Hamiltonian guar- 
antees the existence of unitary time evolution operator exp{—iHt). Due to the 
existence of the black hole, Dirac fields asymptotically approach either horizon 
or infinity. Then we expect, for some data the Dirac fields should approach free 
fields with a free Dirac Hamiltonian Hqo at infinity in distant past and future. 
This means that the following limit should exist 

lim II e"^-^*$ - e-*-^°°*$oo ||= 0. (19) 

i— »±oo 

However, it turns out that the above limit does not exist at infinity for massive 
Dirac fields without certain modification. The possible remedy is to add a new 
phase shift 5{t) to the free wave so that the following modified limit does exist 
at infinity 

lim II e-*-^*$ - e-*^-*e^'^(*)$oo ||= 0. (20) 

To understand this situation better, we like to provide a physical picture before 
we get into tedious mathematical proof. 

3. Physical Picture 

Suppose a black hole is located at O, and at t = a classical Dirac spinor field 
is concentrated around S outside black hole. Denote the distance OS = |xo|. For 
simplicity, let us first assume this free field can be described by a peaked wave 
with a peak energy E. This peaked wave propagates from S at t = to near the 
surface of a sphere with a radius | vt| at a later time t, where v = p/-B is the peak 
velocity. After a sufficiently large time, the closest distance between the wave and 
the black- hole origin is |vt| — |xo|, and the longest distance is |vt| + |xo|. Then the 
propagated wave e~* °°*<l>o(0,x) at a later time t is almost concentrated inside a 
certain domain 

|vt| - |xo| < |x| < |vt| + |xo|. (21) 

Usually, the classical wave operators from (19) are defined by 

VF±$o = s- hm e*^*e"^^-*$o- (22) 

£—* lb CXD 

To prove their existence, by Cook's method it suffices to show 

N{t) =11 {H - H^)e-'"^'^o ||G L\Rt). (23) 

But this is frustrated by a long-range mass term since from (16) 

CMm 
H-H^^-{as0l) + 0{r-^). (24) 



By the domain of the propagated wave (21), it suggests a lower bound for 

Nit)>-^^\\^o\\^Oi\tn. (25) 

\Vt\ + |Xo| 

It turns out that the fower bound of the norm N(t) is of the order 0(|i|~^), which 
is not integrable in L^{Rt) . So the condition (23) is not met in the presence of 
fong range mass term, and the classical wave operators defined by (22) do not 
exist at all. 

If we can somehow offset this long range term by implementing a phase shift 
from the free dynamics, we might construct some wave operators that exist at 
infinity. The modified wave operators from (20) are defined by 

Wt^o = s- lim e^-^'e-^^-V-^^^^^o- (26) 

t^±oo 

It then suffices to show 

Nsit) =11 [H-H^ + ,5'(t)]e-^^-*e*'^W<I>o ||G L\Rt). (27) 

Suppose the phase shift can be adjusted so that 

H-H^ + 6'{t) -^ {-}--)+ 0{t-^). (28) 

|vt| r 

By (21), it suggests an upper bound for 

N^{t) ^11 (Mj_^)e"^-f^-*e^'^W<5o II +0(1*1"^) 

^ I ,m '^°' I n II ^0 II +0{\t\'^) - 0(|t|-2), (29) 

|vt|(|vt| - |xo|) 

which belongs to L^{Rt). Then the condition (27) is satisfied, and the modified 
wave operators defined by (26) do exist at infinity. 

By considering the propagation of a peaked wave, we have shown roughly the 
nonexistence of the wave operators (22) and the existence of the modified wave 
operators (26) at infinity. This gives us a sense what is going on in this kind of 
scattering problems. With this picture in mind, we now go on to provide a strict 
mathematical proof. 

4. Mathematical Proof 

Generally, the classical wave operators are defined by 

W^^oo = s- lim e*-^*Je-*-^°°*$oo, (30) 



where J is a. smooth bounded identifying operator from Tioo to 7i to be specified 
later. To prove the existence of wave operators, by Cook's method [17] it suffices 
to show 

N{t) =11 {HJ- jH^)e-'^-''^oo \\h^ L\Rt). (31) 



As discussed previously, this is frustrated by a long-range mass term 

(^ Mm 

Hj-jH^^-j{a^(^I)^^^^^ + 0{r-^), (32) 



and the condition (31) can never be satisfied. We need to develop a scattering 
mechanism to get away from the Coulomb mass term by introducing a phase shift 
from free dynamics. 

Here the modified wave operators are defined by 

W^^^ = s- lim e'^'je-'^^'e'^'-'^^oo- (33) 

i— +±oo 

To prove their existence, it suffices to show 

Ns{t) =11 [HJ- JH^ + ^(t)]e-^^-*e'^W<I>oo Ht^G L\Rt), (34) 

where 

Sit) = je-^H^'^e'"-\ (35) 

under a condition 

m),6is)]=0, (36) 

for any real t and s. If we choose 

-F/ N / ^N GMm , „, 

<5(t) = :7(a3 0/)— — , (37) 

to offset the long range mass term in (31), then the phase shift 

6{t) = f e^^-*'(a3 ® /)^^e-^-*'dt', (38) 

does not satisfy (36), and is not self-adjoint since [H^,a^ /] 7^ 0. We need to 
find a phase shift that commutes with .ffoo- 

First we need to prove the self-adjointness of the Hamiltonian (18). Let us 
consider a Hilbert space 

Ur, = L^{R X 5^ C^] dr^ded<j)) = L^{R; dr.,) L'^{S^; dOdc/)) C^. (39) 

Then Hr, is a symmetric operator on Cq^{R x S"^; C^) C Tir,- We want to show 
it is essentially self-adjoint (e.s.a.) on this domain along the line of [4]. 



Theorem 1: Hr, is e.s.a. on C^{R x S'^;C'^) C Hr,. 
Proof: The Hamiltonian Hr, in (18) can also be written as 

Hr, =ai(g) H^^ + (0-3 (g) /)m^^/^ 

d fi/2 

H', =-i^—ai + ^K. (40) 

or^ r 

Since ai is symmetric on C^ and (o"3 I)m^^''^ is bounded, by the Kato-Relhch 
theorem [17] it suffices to show H'^.^ is e.s.a. on Cq°{R x S'^;C'^) in H'^.^ = 
L?'{R X S'^,C'^). The operator K of (17) is an angular operator on L'^{S'^;C'^), 
and is e.s.a. on C^ {S"^ ; C"^) , by a similar argument in [18]. Also K has point 
spectrum of finite multiplicity: KZ^ = kZk,k = 0, ±1, ±2, ... [14] [2], where 
Zk are C°° functions. Let Yfc be the eigenspace of K with eigenvalues ±k: 



Yk = {f & L {S , C ) : Kip = ib/c(/9}. Since ctiK = —Kai, we know ciZ^ = Z_ 
and hence o"i is symmetric in Yfc. Let 7ik = L'^{R) ® ^k- We define 



k 



d £1/2 

Hk = -i^r- ®cri + kl. (41) 

or^ r 

In [19], there is a proof of essential self-adjointness for the Dirac operator on 
Minkowski space by Fourier transforming it into a multiplication operator in 
momentum space. This can be applied similarly to show our radial operator —idr, 
is e.s.a. on C^{R). Meanwhile ui is symmetric and ^^/^r"^ (gi kl is bounded, by 
the Kato-Rellich theorem we know Hj. is self-adjoint. Let H'^ = ®kT^k- Then 
H'j.^ = ®kHk is self-adjoint on 

i?«) = {/ : A G D{Hk),Y. II Hkfk f< 00}. (42) 

k 

This extends H'^.^ on the domain V consisting of finite combinations of / (8> Yfc with 
/ G Cq°{R). Now Cq°{R) is a core for —idr^ and hence a core for Hk- It follows 
that I? is a core for H[. , namely, H[. is e.s.a. on V. The symmetric extension of 
H'^^ to the larger domain Cq^{R x S'^; C^) is also e.s.a. with the same closure [4]. 
This completes our proof. □ 

Now making a new coordinate transformation 

^=r^/^(r■) = (l-^-^/^ (43) 

dr r 

we can express the new variable vq by an integral 

ro(r) = r r^^\r')dr'. (44) 

J2GM 



Implementing a substitution of variable 

/ o^^.cosh^x ^, ,, ^ 2GM 1 

r' = 2GM—-^, e(r') = l — = ^-^, (45) 

Sinn a; r cosh x 

we integrate by parts 

C^l\r')dr' = 2GM /(coslia;)d(''°'^' ^^ 



sinh X ' 
, cosh X f cosh X 



„^,^rCOsh X /■ cosh a; ,, , ,, 

= 2GM[ 2 / 2 (i(coshx)] 

sinh X J cosh x — 1 

„^,^rCOSh^X , 1 /", 1 1 N ,/ , X 

2GM 2 coshx / ( )d{cos\ix) 

sinh X 2 J cosh x — 1 cosh x + 1 



coshx 1 coshx + 1 
sinh^ X 2 cosh x - 1 ' 
By (45), we write 

ro(r) = re^/2 + GM log \'^_%r (47) 

When r = 2GM, we know ^^ = and ro(2GM) = 0. When r -^ oo, we have an 
asymptotic form 

ro|,^oo = r + GMlog(r) + 0(l). (48) 

This asymptotic form can also be derived by integrating (43) 

^1™ = r'/\-.oo - 1 + ^ + 0{r-'). (49) 

dr r 

By taking dro/dr^, = ^^'^, we rewrite Hamiltonian (18) as 

Hro = ^1 kl(-^e'/'T^) + —K] + (dg /)mCV2^ (59) 

with a Hilbert space 

Hro = ^^((0, 00) X S^; C^; C^^^drodedcp). (51) 

Lemma 1: Hr^^ is e.s.a. on Cq°{{0,oo) x S'^;C*^) C Tiro- 

Proof: From the coordinate transformation r* -^ rg, we define a unitary 
transformation Ui : Hr, -^ Ti-ro t>y 

{Ui'i>){ro,e,(^) = <i>{r,iro),e,^), (52) 

This is unitary since dr^/dr^ = ^^'"^ and so 

f{Ui<P)^{Ui^)r^/^droded4>= j ^^^dr^dOdcj). (53) 



In general, unitary transformations preserve self-adjointness. We conclude H^q 
is e.s.a. on smooth functions with compact support in Htq since H^, is e.s.a. on 
smooth functions with compact support in Hr^ by Theorem 1, and since Ui maps 
C^ to C^. □ 

The free Hamiltonian can be written asymptotically tq ^ oo as r ^ oo 

d 1 

H^ = ai® [ai(-i-— ) + —K] + (^3 ® I)m, (54) 

dro ro 

with a Hilbert space 

Hoc = ^^((0, cx)) X S^; C*; dr^dedcf)). (55) 

Proposition: Hoo is e.s.a. on Co°((0, 00) x S'^;C^) C Hoc- 
Proof: It is a special case of Lemma 1 when M = 0. □ 
From (50) and (54) we have 

HrgJ— jHao = J{i — l)Hoo 

+J{ai K)e/'{- --)- (cTi aiM'/^^). (56) 

r ro dro 

Here the identifying operator J is defined for some 6 > a > 

je C°°iR), < Jiro) < 1, (57a) 

J{ro) = for a > ro > 0; J^(ro) = 1 for tq > b. (576) 

From (48), the second term in (56) has an asymptotic order 

e/\- --) = ^^^ + 0{r--) = 0{r-^^n = 0{rf^% 



ro r^ 



(58) 



for log(r) < C^r^ with e > 0. The third term in (56) vanishes as r ^ cx) since 
dj/dro = for ro > b. In the first term of (56), we use (58) to estimate 

^1/2 _ 1 = _™ + o{r-^) = -^ + 0(ro-2+^). (59) 

r ro 



To offset this term, let us choose 



^ii) = ^J^^--- (60) 



where |v| = \J—/\/\J—/\ + m,^ and A is the Laplacian in spherical coordinates. 
Then from (35) we find a logarithmic phase shift 

5{t) = ^H^\og{t), (61) 



10 



where 

/ s , , , dlos,(t) 1 , , 

log(t)=s<7n(t)log|t|, —^ = - (62) 

By changing Hilbert space, we now work directly on the new coordinate tq to 
present our main theorem. 

Theorem 2: For the classical scattering of massive Dirac fields off a black 
hole characterized by a Hamiltonian: 

ovq r 



with C = 1 - 2GM/r and 



9 15 
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the modified wave operators defined by 

Wt^oo = s- lim e*^'-o*Je-^-f^°°*e*''W$oo, (65) 

do exist in T^^o by picking a time-dependent logarithmic phase shift 

5{t) = ^H^log{t). (66) 

Proof: Since the time-dependent operator in (65) is bounded, it suffices to 
show the existence of the limit in (65) for <I>oo in a dense set Di. Di is defined 
so that in Cartesian coordinates, the Fourier transform has compact support and 
vanishes in a neighborhood of the origin. Its exact definition will be given in 
Lemma 2. Then by e/3 arguments, we know the existence of (65) for ^^o in the 
whole Hilbert space. By Cook's method, the existence problem is reduced to 
show 

Nsit) =11 [Hr,J- ^i^oo +^(t)]e-'^-*e^^W<I>oo ||h.„G L\Rt). (67) 

To obtain this, we have to make sure that Je~ °°*e* '*^<I>oo is in the domain of 
Hrty We leave this problem to Lemma 2. 

With the logarithmic phase shift (66), we may estimate the norm in (67) by 
considering the asymptotic property of each term in (56) 

Nsit) <\\ ^(^ - ^)$t 11^^^ + II JZiro)'^t 11^^^, (68) 

where Z{ro) is bounded and Z{ro) = 0{rQ ^'^) as ro -^ 00, and also 

Mt,ro,0,(p) = GMi7ooe-'^°°*e^^W$oo(ro,^,0). (69) 

11 



To change Hilbert space, we estimate 

< CiJlJFfdrodedcpy/^ = C\\JF ||^^< C \\ F ||^^, (70) 

where we have used the following inequalities 

1^2^-1/2 < |^2^2 < ^2^ (7^) 

with an upper bound C^ = ^max = ^^^'^{ro\min)- So (68) becomes 

Ns{t) < C II ^(^ - !)<&, 11^^ +C II Z{r,)^t \\n^ . (72) 

Transforming back to Cartesian coordinates 

xi = ro sin cos 0, X2 = tq sin ^ sin 0, X3 = rocos^, (73) 

we get a matrix \ul> defined by dxi = uldyj;i,j = 1, 2, 3 where yi = ro, y2 = 
and y3 = (/). Denote the determinant g = det{g'^'^) = det^g"^^). To preserve the 
invariance of space-time interval, it is easy to get (7(x) = g{y)/\det{u)\'^ where 
det(ti) = r^(sin0). Then we get a transformation for half density 

In Minkowski space-time, $oo = (— 5')^ ^^oo; and (—5)^'^ transforms by (74), so 
we can define a coordinate transformation operator T(u) as 

(r(n)$oo)(x) = (^-|^)(ro(x),^(x),</>(x)). (75) 

We rewrite the free Hamiltonian from (54): 

i^oo(ro,^,0) = il'h'^ + -^ + ^^tS;) + ^°"^' (76) 

otq ro off ro sm ocp 

and make a transformation T[u) by substituting coordinates 

THoo{ro, 0, 4>)T^^ = ij^Kj^ sin^cosc/) + 7^ cos0cos(/) — 7^ sin 

d 

+ (7 sin0sin(/> + 7 cos ^ sin </> + 7 cos (/>) — — 

0x2 
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d 

dxi 



d 

+(7^ cos 6* -7^ sin 6*)- — 1 + 7°m = H'. 
0x3 



To transform (77) into 



,0. .1 ^ 



d 



d 



//oo(x) = i7"(7'T^ +7'i^ + 7'i^) +7""^ 



dxi 
we find a spinor rotation defined by 



dxo 



dx^ 



(77) 



(78) 



sm & cos (p sm t> sm (p cos t> 

cos cos (^ cos 9 sin — sin ( 

— sin 6 cos d> 



(79) 



Here R is orthogonal for R ^ = R and det(i?) = 1. The spinor representation 
S{R) of such a rotation is unitary S~^ = S^ [13]. We may now write 

SH'S-^ = //oo(x). (80) 

Let us rephrase the above problem by defining a dreibein field 

JV^^I = diag{l,r-^,ir sin ey^). (81) 

Then (76) can be written as 



Hoo{ro,(^,i 



nVv^'9, +7V, 



(82) 



where {dj} = (drt) , de , d^) . Making a coordinate transformation dj = (u )j(9^, 
where {i9[,} = {dxi,dx2idx3), we get the equation (77): 



(83) 



Define a new dreibein field (V)^ = V-' {u )^ in Cartesian coordinates, and a 
rotation matrix R^^ = I {V')^ > where the elements of R is given by (78). Under 
a spinor rotation S{R)'y'^S~^{R) = Rj'y^ we arrive at the equation (80): 

SH'S-^ = i7°Vi?j(^"^)f ^fc + 7°"i = ^7°V^; + 7°"^ = ^oo(x). (84) 

Finally we define a unitary transformation as a product of the coordinate 
transformation and the spinor rotation 



U = T{u)S{R), 

mocK^) = ( |ft^i^r/2 )(^o(x),^(x),0(x)). 



(85a) 
(856) 
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This operator preserves the norm squared 

S{R)<^ 



<^oo{ro,e,<P)\^droded<p. (86) 



And U is unitary from Hoo to 

n'^ = L\R^;C'';d^^). (87) 

Making this unitary transformation in (72) and defining ^t = U^t, we get 

Ns{t) < C II J(|x|)(^ - -^)^t \\n' +C II ^(|x|)^i 11^' , (88) 

with 

^i(t,x) = GMi7ooe-'^°°*e'^(*)^oo(x), (89) 



where |v| = y/—A/\/—A + tti^ and Hoo = —ia • V + Pm in Cartesian coordi- 
nates. Note here the norms in (88) are taken on L'^{R^;C'^), while those in (72) 
onL'^{R+ xS^;C^). 

Since Hoo has a symmetric energy spectrum, any spinor can be expressed by 
a superposition of the positive and negative energy components: 

*oo(x) = ^i(x) + M/- (x), (90) 

where 

^±(x) = |u±(p)e±^P-d3p, (91) 

with spinors satisfying 

{a-p + pm)u±{p) = ±EpU±{p), (92) 

Then the transformed wave (89) sphts into two parts as well 

*i(t,x) = ^+(t,x) + ^,-(t,x), (93) 

where 

*,±(t,x) = ±y'GMi^pt/±(p)n±(p)e±^P-^d3p, (94) 

with two time-dependent unitary operators 

Utip) = expi^iEpt) exp[±iEp^logit)]. (95) 
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By the Minkowski inequality, we have 
Nsit) < C II ^(|x|)(^ - ^)*+ 11^, +C II J{\^\)i-^ - j^)^t \\h' 

\\'t\ |X| "-°c Ivtl |X| "-oo 

+0 11^(1x1)^^+11^, +C II Z(|x|)*,- 11^, . (96) 

Since the unitary operators U^ (p) are of the type discussed in [3] and each 
component of (94) is a smooth solution of the Klein-Gordon equation with the 
Fourier transform in Co°(i?^\ {0}), all four terms on the right-hand side of (96) 
can be proved to belong to L^{Rt) by the general method provided in [3]. 

Since clearly .Z^(|x|) = 0(|x|~^"'''^) G L'^{R^), the last two terms on the right- 
hand side of (96) are of the order 

II Z(|x|)^± 11^, <|| Z(|x|) Ibll ^t lloo< 0(ir'/'+^), (97) 

which belong to L}{Rt) by the theorem 1(a) in [3]. By the definition of the 

identifying operator J (57), the first two terms on the right-hand side of (96) are 

of the order 

II ^(|x|)(^ - ^)^t \\^, < 0(|t|-2+^), (98) 

Ivil |x| "-00 

which also belong to L^{Rt) by the theorem 1(b) in [3]. 

With these results, we arrive at the condition (67) Ns{t) € L^{Rt), and thus 
prove the existence of the modified wave operators (65) with a time-dependent 
logarithmic phase shift (66). As in a Corollary in [3], it is easy to show the 
modified wave operators W^ are independent of the smooth bounded identifying 
operator J. □ 

Lemma 2: Define a dense domain 

Di = {<peHoo. {U<t>r e C^{R^\{0}; C*) c s], (99) 

where f^ is a unitary transformation defined by (84), {U(j))^ is the Fourier trans- 
form of U(j), and S is the infinitely differentiable Schwarz space. Let <J>oo G Di- 
Then Je~*^°°*e*'^(*)<^oo is in the domain of Hr^y 

Proof: Since H^ does not change Di C TYoo, with (66) we know 

if = exp[-iH^{t - ^ log(i))]'^oo G Di. (100) 

By the previous Lemma 1, i/j-o is e.s.a. on Co°((0, oo) x S'^;C'^) C D{Hr^^. To 
show 7/> = Jip G D{Hrf)), it suffices to show there exist ipn G Cq° so that Hr^ipn 
converges when ipn ^ ip- Define a series of characteristic functions Xn for some 
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small ai and bi: < ai < bi < a, where a is the same as the one used in defining 
the characteristic function J (57): 

XnGC°°(ii), 0<Xn(ro)<l, (101a) 

Xnifo) = for < ro < ai or ro > n + ai, (1016) 

Xniro) = 1 for hi <ro < n. (101c) 

The characteristic function Xn cuts off wave function outside {ai,n + oi). If we 

define ipn = Xn^P G C^, then we know ipn -^ Xooip = ip and 

Hrolpn = 0-1 (g) 0-l(-^^^^^)(^^V' + Xr; 



oro orQ 



£l/2 

^ ki ® [(Ti{-ii^'^^^) + ^K] + (as /)me^/V, (102) 

aro r 

where we have used the result that {dxn/dTo)ip -^ pointwise. Also dxn/dro < 

/(ro) G L^, by the dominated convergence theorem, we know || {dxn/drQ)ilj \\^ 0. 

By (102), Hrf^Tpn converges to -ffroV'; we thus conclude ip £ D{Hrf^). □ 

5. Remarks 

Technically, the choice of a new coordinate ro rather than r^, is of critical help to 
find a phase shift that commutes with the free Hamiltonian. The time-dependent 
logarithmic phase shift we have derived is a matrix operator due to the mix-up 
of positive and negative energy spectrums. As a result, there is a sign difference 
between the phase shifts for pure positive and negative energy waves. 

With the results of this paper about massive spinor fields and those of [3] 
about massive scalar fields, we can now conclude that for the scattering of some 
massive fields by black hole it is possible to find a logarithmic phase shift from 
the free dynamics to offset the long range mass term at infinity. 

If we also consider a local interaction V{x), we need to replace the free bare 
mass m by a locally dressed mass m + V{x) [2] in the Dirac equation. The proof 
of the existence of modified wave operators at infinity in this paper should still 
be true as long as V{x) -^ reasonably fast at infinity. 

Although we choose a simple spherical black hole as our manifold, we can also 
apply our theory to any spherical stars with a Schwarzschild metric outside. 

In addition to scattering, there are still many interesting problems that remain 
to be clarified in the relativistic theory of quantum Dirac fields on black hole such 
as the Hawking effect. 
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